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Abstract 

A universal description is proposed for generic viscoelastic systems with a single relaxation 
time. Foliation preserving diffeomorphisms are introduced as an underlying symmetry which 
naturally interpolates between the two extreme limits of elasticity and fluidity. The symmetry 
is found to be powerful enough to determine the dynamics in the first order of strains. 
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1 Introduction 



Viscoelasticity is a notion that unifies solids and fluids (see [T] for example), and is applied to a 
wide range of materials as in rheology [2]. Viscoelastic materials behave as solids during short 
intervals of time, while they do as ordinary viscous fluids at long time scales [3JH]. 

In order to get a concrete image for viscoelastic materials, let us imagine something like a 
chewing gum. We can see many particles bonding in it to each other (see Fig. [I]). When the 
material is stressed, the bonds produce an elastic force and try to make the particles back to 
the original configuration. That is, the system exhibits elasticity during short intervals of time. 
However, if we keep the deformation for a long time, then the bonding structure changes to reduce 
the free energy as in Fig.[H and the shear stress vanishes. We further assume that the material is 
elastic for compressions even for long time scales, so that the bulk stress does not undergo such 
relaxation (see Fig. [2]). Thus the system exhibits fluidity at long intervals of time. 

shear 
deformation 
>■ 

Figure 1: Shear deformation of viscoelastic material. If we keep the deformation for a long time, 
the bonds are reconnected and then the stress due to the deformation vanishes. 






Figure 2: Compression of viscoelastic material. The material we consider here is elastic for 
compressions even for long time scales, and thus the bulk stress does not undergo relaxation. 



In this paper, we propose a framework in which the characteristics of elasticity and fluidity 

can be both dealt with on equal footings. We introduce two key ingredients. One is a new variable 

(to be called the intrinsic metric) which represents the bonding structures. The other is the gauge 

symmetry of foliation preserving diffeomorphisms (FPDs)[j] which we find to interpolate between 

the two extreme limits of elasticity and fluidity in a natural manner. We show that requiring 

the invariance under FPDs uniquely determines the dynamics of viscoelastic systems in the first 

1 An interesting application to quantum gravity was recently found by Hofava [5]. For mathematical details of 
FPDs, see also [6]. 
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order of strains. 



2 Geometrical setup 

For an elastic material, a sufficiently small region around any point can always be regarded as 
being deformed from the configuration with no strains. We call the shape before the deformation 
the natural shape (see Fig.[3|). This is a straightforward generalization of the notion of the "natural 
length" for a spring or a rubber string. For a viscoelastic material, we assume that it exhibits 
the elasticity during short time intervals, so that we can define the natural shape around a given 
point at each moment. Since the natural shape is constant in time for elastic materials, its time 
dependence represents the plasticity (i.e. non-elasticity) of the material. In order to describe the 
natural shape quantitatively, we introduce a new dynamical variable to be called the intrinsic 
metric, and discuss its properties in this section. 



natural shape 



Figure 3: A sufficiently small region around any point can always be regarded as being deformed 
from the natural shape. 

At each time t, we introduce coordinates £ = (£ a ) (a = 1,2,3) arbitrarily on the material^] 
We then define the intrinsic metric of the material (denoted by g a b(t,,t)) such that the distance 
between two points at £ and £ + d£ at fixed time t is given by the distance in the natural shape 
(see Fig. ED : 

ds 2 = g a b(£,t) d^ a d^ b = (length in the natural shape) 2 . (2-1) 

We let X*(£,t) be the spatial Cartesian coordinates of the point with £ at t . We emphasize 
that the intrinsic metric g a b(^,t) is independent of X l {^,t), and can differ from the induced 
I 



metric 



h ab ti,t) = d a X\Z,t)d b X\Z,t), (2.2) 



2 Here we assume that the coordinates £ move smoothly on the material as the time t varies. 

3 For a function F(f,i) we write F = dF/dt and d a F = dF/d^ a . 
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natural shape 

Figure 4: The intrinsic metric is denned at each time t as the distance in the natural shape, ds 2 
Qab{i,t) d£ a d(; b . The real distance is measured with the induced metric, ds 2 ^ = h a b(£,t) d£ a d£ b . 



which measures the distance in the real three-dimensional space (see Fig.0]). Note that their 
discrepancies 

£o&(£,*) = 5 (M£, *)-&&(£,*)) (2.3) 

represent the strain tensor of the system. In order to see this, let us consider an elastic material. 
We take the coordinates £ such that they move with the atoms and coincide with their Cartesian 
coordinates in the mechanical equilibrium. When the material is deformed slightly, the atom 
labeled with £ has the coordinates X a (^,t) = £ a + u a (£,t), where u a (^,t) is the displacement 
vector, and we have 

Kb = S ab + d a u b + d b u a + 0(u 2 ) . (2.4) 

Since the intrinsic metric keeps taking the values in the mechanical equilibrium, g ab = 5 ab , we 
have that e ab = (l/2)(h ab — g ab ) ~ (l/2)(d a u b + d u a ). This shows that e ab is certainly the strain 
tensor for elastic systems. In the following discussions, we make the linear approximation with 
respect to £ ab ■ This implies that for the quantities of 0(e ab ) one can raise or lower vector indices 
by either of g ab or h ab . 

We denote by po the mass density in the absence of strains. In this paper, we assume that po is 
a constant independent of £ and t. Then the mass contained in a volume element d 3 £ = d£} dt^ 



is given by poy/ g{£, t) <i 3 £, and the mass density in the real three-dimensional space is given by 

As an example, let us consider a squeeze deformation of a two-dimensional viscoelastic material 
formed by bonding particles as in Fig. [5]. The wavy lines stand for the bonds. The left figure 
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Figure 5: Squeeze with the angle 9 . This can be obtained by combining a shear deformation and 
a compression (plus a rotation). 

represents the material before the deformation. Here we take £ to be Cartesian coordinates, 
and then both of the induced and intrinsic metrics have the same form: h a b = g a b = $ab ■ 
Therefore the mass density p is equal to pg • We assume that £ are attached to the particles, 
so that they comove with the material under the deformation. The middle figure represents the 
material just after the deformation. The induced metric and the mass density change according 
to the deformation, while the intrinsic metric does not. The induced metric is evaluated as 
ds 2 = (dt; 1 ) 2 + (d£, 2 ) 2 — 2d£}d£, 2 cos(7r/2 + 9) by using the cosine formula. The right figure 
represents the material after a sufficiently long time. In the process of relaxation, the induced 
metric and the mass density are preserved, while the intrinsic metric g a b becomes proportional 
to h a b , g a b = /(£, t) h a f, , as will be discussed around (|4.7p . The factor /(£, t) can be determined 
from p through (|2,5|) . 

3 Foliation preserving diffeomorphisms 

Since we have introduced a different coordinate system £ at each time, in order to describe the 
actual motion of a material one needs to specify how the fluid particles move relatively to the 
coordinate system £. This can be realized by introducing a vector field N a (^,t), with which the 
fluid particle located at £ at time t is supposed to move to the position £ a + N a (£,t)5t after the 



4 



time interval 5t (see Fig. [6J. 




Figure 6: Definition of N a 



In fact, by using N a , the time derivative along the path of a fluid particle (the material 
derivative) of a scalar quantity <p(£,t) can be expressed as 

^=<p + N a d a <p. (3.1) 

For example, when measured in the real space, the velocity v l (£,t) and the acceleration a l (£,t) 
of the fluid particle located at £ are given by 

* DX i 



Dt 

Dv { 



X l + N a d a X\ (3.2) 
u < + JV°a o « < . (3-3) 



Dt 

Then we consider the foliation preserving diffeomorphism (FPD): 

t' = t, £'° = £'*(£, t). (3.4) 

We say that a quantity T°"L (£,<) is a covariant tensor, if it transforms as a three-dimensional 
tensor at each time: 

T >a ^',t) = ^---^ b ---T c -- d ..X^t), (3.5) 

where the transition functions d^' a /d^ c and d£ d /d£,' b should be evaluated at time t. For example, 
X 1 , v l and a 1 are covariant scalars, and g a b, h a f, and e a b are rank two covariant tensors. 

On the other hand, N a is not a covariant vector but transforms with an inhomogeneous termQ 



4 This is naturally understood, if one recognizes that N a can be regarded as a gauge field that appears when 
the three-dimensional diffeomorphism is gauged in the time direction. 
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To see this, let us consider a fluid particle in two different coordinate systems, and assume that 
it has the coordinates £ a and £ /a , respectively, at time t. Then, by definition, at time t + St, 
they become £ a + N(£,t)5t and £ /a + N /a (^' ,t)St, respectively. Because these are related by the 
transformation (|3.4p at t + St, we have 

+ iv /a (£', t)st = e' a (e + n(c, t)st, t + st) 

and thus we obtain (|3.6p . 

Note that the time derivative of a covariant tensor T a " fe ... is no longer covariant. We, however, 
can make a covariant tensor by considering the derivative along fluid particles as in (|3.2|) and 
{331: 

^ •/•""/,.. •/•"••/,.. ' £ A - -r'",,.... (3.8) 
where £/v is the Lie derivative with respect to the vector field M = N a d a : 

£ M T a ~ b „, = N c d c T a ~ b ... - d c N a T c ~ b ... + d b N c T a ~ c ... + ■■■ . (3.9) 

For example, from g ab , we can make a covariant tensor K ab to be called the extrinsic curvature, 

Kab = \ iff = \ ^ ab + ^ 9cb + &hNC 9ac + ^ 9c9ab> • (3,10) 
Similarly, we can define the extrinsic curvature which corresponds to h ab , 

K ^ = \^r = \ + d * NC h <* + dbNC h ^ + NC ■ ( 3 - u ) 

Note that 

K^-K ab = ^. (3.12) 
We then introduce the dreibeins e l a (^,t) and their inverses e"(£, t) : 

e^,t) = d a X i ^,t), e?ei = 5Z. (3.13) 
From (|3.2p and (|3.13j) . we obtain 

h ab = eiet, d a el = d b el, ^ = d a v l - (d a N b )et - N b d a e\ , (3.14) 

and thus 

-pf = (3.15) 
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By using the dreibeins, the scalars v l (£,t) and a l (^,t) can be transformed into vector fields as 

a a i i i a a i i i /o 1 r\ 

v = e^v , v a = e a v , a = e^a , a a = e a a . (3.1b) 

We can calculate the acceleration a a by using the Leibniz rule for D/Dt and f)3. 15|) as 

Dv* D(e^) >Dei Dv a . /l,^ 



= v a + (d a N b )v b + N b d b v a - d a (]^ b v h \ . (3.17) 
Furthermore, by using the first of (|3. 14|) and (|3.15p . we can show 

= \ 4 + 4) = \ (V a h) v b + Vf } «.) • (3.18) 

The set of FPDs forms a gauge symmetry group of the system, which can be gauge-fixed 

arbitrarily according to convenience in describing the dynamics of a given system. Two of the 
useful gauge fixings are the followings: 

(A) comoving frame: 

We set N a (£, t) = . In this frame, X*(£, t) describes the motion of the fluid particle attached to 
the coordinate £ , and we have 

yifa t) = JT(£, t) , a l {C, t) = X\£, t) (in the comoving frame). (3.19) 

This frame is useful for describing the dynamics of elastic materials. Note that the acceleration 
field (|3.17p is expressed as 

a a = Va — d a ( — v b v j (in the comoving frame). (3.20) 



2 

The last term gives the inertial force. 
(B) laboratory frame: 

We set X a (£,t) = i a . In this frame, we have X*(£,t) = 0, ef = 5°; and h ab = S ab ■ Then, from 
(EOT) and mngp we obtain 

v a (£, t) = N a (£, t) (in the laboratory frame). (3.21) 

This frame is useful for describing the dynamics of fluid. Note that the acceleration field (|3.17|) 
represents the material derivative of the velocity field 

a a = v a + v b d b v a (in the laboratory frame). (3.22) 
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4 Fundamental equations 



We are now in a position to write down a set of equations which determine the time evolution of 
X l (£,,t), g a b(£,t) and N a (^,t) up to FPDs. The covariance under FPDs is found to be powerful 
enough to uniquely determine the dynamics in the first order of e a b ■ 

We first consider the limiting case of elasticity. Due to our definition of the intrinsic metric, 
9ab(€it) should be constant in time for elastic materials in the comoving frame (N a = 0). The 
FPD-covariant expression for this statement is that the extrinsic curvature K a b of (|3.1U|) vanishes: 

K a b = (elastic limit) . (4.1) 

This implies that the non- vanishing K a b represents the genuinely plastic deformations. 

In order to make further discussions, we need to separate the trace part from K a b because it 
vanishes for any materials due to the mass conservation of the system: 

K(U)^9 ab (^t)K ab (C,t) = 0. (4.2) 

In fact, contracting (|3.10j) with g ab , we obtain 

K = \ g ab {g ab + d a N c g cb + d b N c g ac + N c d c g ab ) 

= -L((Vgy + d a (^g-N a )) 

= — {pint + d a ( Pint N a )), (4.3) 

Pint 



where /9mt(£> t) = poy/g{£, t) is the mass density with respect to the intrinsic metric. This indicates 
that the vanishing of K is equivalent to the mass conservation, which can be easily seen in the 
comoving frame where N a = and p m t = 0. 

The traceless part of K ab on the other hand describes the rate of the shear deformation of the 
intrinsic metric, and thus is expected to be proportional to the traceless part of the strain tensor 



Kab(Z,t) = -e ab (t,t), (4.4) 
r 



where r is the relaxation time, and 



Kab = K ab ~ ^Kg ab , (4.5) 

Sab = £ab ~ g {g Cd £cd) gab ■ (4.6) 

We call (|4.4|) the rheology equation hereafter. After a time interval much longer than r, £ a b(£, t) 
vanishes and thus g a b{£,,t) becomes proportional to h a b(£,,t): 

9ab(£, t) = /(£, t) h ab {£,, t) (/(£, t) : a scalar function) , (4.7) 
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where /(£, t) is determined by the mass conservation as in the example in Sec. [2J 

The equation (|4.4p is the simplest and is expected to be universal. It is consistent with (|4.1D 
because the elastic limit corresponds to r = oo . The fluid limit is also realized correctly by taking 
r = , in which we have i a b = . This implies that g a b is proportional to h a b as in (|4,7p . The 
only remaining degree of freedom of g a b becomes the density of the material, which means that 
the system corresponds to an isotropic fluid. 

Note that (|4.2p and (|4.4p indicate that all the components of K a b are of the order e. Then 
from (|3.12p we find that so are those of and we have 

K^-K ab = ^ + 0(e 2 ), (4.8) 
K 00 - K = + 0{e 2 ) . (4.9) 

The dynamics of X l (£,t) should be expressed as Euler's equation which is written with the 
stress tensor T a b as 

pa a = -V^ b T ba , (4.10) 

where is the covariant derivative with respect to h ab ■ We show that the leading form of the 
stress tensor in the derivative expansion can be determined by the following requirements: 

• T ab is symmetric and covariant under FPDs. 

• T a b is linear in the strain e a b and the spatial derivative of the velocity, Va^v&. 

The above requirements imply that T a b is a linear combination of the irreducible components 
of E ab and K$> = {l/2)(Va h) v b + V ( b h) v a ) (see §MY- 

Tab = -2ne[ h b ] - -eW h ab - 2 1 K^ - QK^h ab , (4.11) 

where 

£ (V = h ab e ab , (4.12) 

e a h b ] ^e ab -\e^h ah , (4.13) 
K (h) _ h ob K W ^ ( 4-14 j 

K^^K^-\ K ^hab. (4.15) 

From the discussions around (|4.8p . we find that all the terms in (14. lip are of the order of e. 
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In order to see the meaning of the coefficients in (|4.1ip . we consider two extreme limits of 
elasticity and fluidity. We first consider the elastic limit. Since K ab = in the elastic limit (see 
(gZED ), the formulas P~8|) and P~9j) lead to 

(h) 



K^ = ^(l + 0(s)), (4.16) 
K= — (l + 0(e)), (4.17) 



and thus we have 

1 Di^ De^ 

T ab ~-2fi ef b > - - e^h ab - 2 7 -^ - C-^afi (elastic limit) . (4.18) 



This indicates that the parameters fi and 1/k are the shear and bulk moduli, respectively. The 
last two terms express frictions. 

On the other hand, the fluid limit is realized by considering the case where the time scale T 
of the variation of the shear strain e ab , 

D 1 , . 

j^e ab ~-e ab , (4.19) 

is much longer than the relaxation time r, T ^> r. We then can show that 

&* = K% ) (l + o(£)) , (4.20) 



because the following holds due to f|4.8j) : 



7-) 7-)~ 

-K ab = ^{l + 0{e)) = -^f (1 + 0{e)) ~ - e ab = -K ab . 



Since e^j = e ab (1 + 0(e)) = r K ab (l + O(e)) , we can rewrite the stress tensor ()4.1ip into 

T ab ~ -27? k£ } - C ^ (h) /iafe - - h ab (fluid limit) , (4.22) 

where 

77 = 7 + /ir. (4.23) 

By using (|3.18p . each term in (|4.22p can be interpreted in terms of fluid mechanics if we take the 
laboratory frame (X a (^,t) = £ a , v a = N a ): 



Kb = <U , (4.24) 
2 (daVb + d b v a ) - - 

R (h) = h ab K (h) = QcVc {4 26) 



Kab = 7i ( d ^b + d b v a ) - t d cVc 5 ab , (4.25) 
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which indicates that r\ and £ represent the shear and bulk viscosities of the fluid, respectively!^ 
The third term in (|4.22p can be interpreted as the pressure: 

p = --e {h) . (4.27) 

K 

To see this, we notice that eW (~ (y/h - y/g)/y/g) measures the deviation of the real volume 
element from that of the natural shaped and that the pressure vanishes in the natural shape. 
Therefore p must be proportional to — We then find that re agrees with the coefficient of 
compression because re = pdp~ l /dp\ (p = po^/g — Po(l — e( h ') and p = —K~ 1 e^). 

The set of equations (|4.2p , (j4.4j) , (|4.10|) , (|4.1ip and (|4.23|> are the fundamental equations which 
govern the dynamics of a given viscoelastic material. We can see that macroscopic properties of 
such materials are characterized only by six parameters (po, r, 77, £, re, 7) from which another 
parameter p is obtained via (|4.23ft . 



5 Conclusion and Discussion 

We have proposed a set of fundamental equations to describe generic viscoelastic systems in a 
unified way. It is expressed as a world volume theory with the target space coordinates X l (£, t) and 
the intrinsic metric <? a b(£> t). FPDs play an important role in interpolating the two extreme limits 
of elasticity and fluidity. We have shown that the covariance under FPDs uniquely determines 
the form of the equations in the first order of strains. We thus conclude that the set of equations 
gives a universal description of viscoelastic systems. 

Here we have considered viscoelastic systems with a single relaxation time. It would be 
interesting to consider a generalization to the case with more than one relaxation time in order 
to describe more realistic materials realized in laboratories. 
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